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In this paper, we demonstrate that optimal control algorithms can be used to speed up the imple- 
mentation of modules of quantum algorithms or quantum simulations in networks of coupled qubits. 
The gain is most prominent in realistic cases, where the qubits are not all mutually coupled. Thus 
the shortest times obtained depend on the coupling topology as well as on the characteristic ratio 
of the time scales for local controls vs non-local (i.e. coupling) evolutions in the specific experimen- 
tal setting. Relating these minimal times to the number of qubits gives the tightest known upper 
bounds to the actual time complexity of the quantum modules. As will be shown, time complexity 
is a more realistic measure of the experimental cost than the usual gate complexity. 

In the limit of fast local controls (as e.g. in NMR), time-optimised realisations are shown for 
the quantum Fourier transform (QFT) and the multiply controlled NOT-gate (c" _1 not) in various 
coupling topologies of n qubits. The speed-ups are substantial: in a chain of six qubits the quantum 
Fourier transform so far obtained by optimal control is more than eight times faster than the 
standard decomposition into controlled phase, Hadamard and SWAP gates, while the C n_1 NOT-gate 
for completely coupled network of six qubits is nearly seven times faster. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Yz, 03.67.Pp; 82.56.-b, 82.56.Jn, 82.56.Dj, 82.56.Fk 



I. INTRODUCTION 

A key motivation for using experimentally controllable 
quantum systems to perform computational tasks or to 
simulate the behaviour of other quantum systems Q, Q 
roots in reducing the complexity of the problem when 
going from a classical setting to a quantum setting. The 
most prominent pioneering example being Shor's quan- 
tum algorithm of prime factorisation 0, ^J. While classi- 
cal prime factorisation algorithms are of non-polynomial 
complexity NP |5j, Shor's quantum algorithm brings it 
down into the class of polynomial complexity P. Another 
celebrated example is Grover's quantum search algorithm 
0, 0] , which allows for searching in an unstructured data 
base of n qubits with N — 2™ items in 0(y/N) quantum 
steps instead of O(N) classical ones. 

As a matter of fact, many quantum algorithms can be 
subsumised as solving hidden subgroup problems in an ef- 
ficient way In the abelian case, the speed-up hinges 
on the quantum Fourier transform (QFT): while the net- 
work complexity of the fast Fourier transform (FFT) for 
n classical bits is of the order 0(n2 n ) HH3, the QFT for 
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n qubits shows a complexity of order 0(n 2 ). For imple- 
menting a quantum algorithm or a quantum simulation 
in an experimental setup, it is customary to break it into 
universal elementary quantum gates [llj . Common sets 
comprise e.g. (i) local operations such as the Hadamard 
gate, the phase gate and (ii) the entangling operations 
CNOT, controlled-phase gates, -y/swAP, i swap as well 
as (iii) the SWAP operation. The number of elementary 
gates required for implementing a quantum module then 
gives the network or gate complexity. 

However, gate complexity often translates into too 
coarse an estimate for the actual time required to imple- 
ment a quantum module (see e.g. |!2lll3l .lT3|). in partic- 
ular, if the time scales of a specific experimental setting 
have to be matched. Instead, effort has been taken to give 
upper bounds on the actual time complexity M 

, which is 

a demanding goal from the algebraic point of view. With 
the time required for implementing a module in a specific 
experimental setting as the most realistic measure of cost, 
here we use methods of optimal control theory to find the 
minimum time by trying to solve the time-optimal con- 
trol problem. The solution is hard to come by in general, 
so here we resort to numerical algorithms. The shortest 
times obtained depend on the coupling topology as well 
as on the characteristic ratio of the time scales for local 
controls vs non-local (i.e. coupling) evolutions and thus 
embrace the specific experimental setting. Relating these 
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minimal times to the number of qubits gives the tightest 
known upper bounds to the actual time complexity of 
the quantum modules in a realistic experimental setup. 

Moreover, as will be discussed, in the generic case there 
is no simple one-to-one relation between time complexity 
and network complexity, because of different time-scales 
between local and non-local controls, different coupling 
topologies allowing for different degrees of parallelisation, 
and different types of coupling interactions matching dif- 
ferent sets of elementary gates. 

Thus here we leave the usual approach of decompos- 
ing gates into sets of discrete universal building blocks. 
Instead, the scope is to exploit the differenti al g eome- 
try of the unitary group for optimisation [lH 1 1 7t | when 
using the power of quantum control in order to obtain 
constructive bounds to minimal time both as close to 
the experimental setting and as tight as possible. In the 
limit of zero cost for the fast local controls (as in NMR) 
compared to the slow coupling interactions, we give de- 
compositions for the QFT and the multiply-controlled 
NOT-gate C n_1 NOT that are dramatically faster than the 
fastest decompositions into standard gates known so far. 

The paper is organised as follows: the first focus is on 
the fact that for time-optimal decompositions of a de- 
sired unitary gate into a sequence of evolutions of exper- 
imentally available controls the global phase may play 
a role. This is the case when e.g. there are different 
time scales for local versus non-local controls. However, 
global phases can readily be absorbed by shifting gra- 
dient flows from unitary to projective unitary groups. 
Then numerical time-optimal control provides the cur- 
rently best upper bounds to the actual time complexity 
of quantum modules like the QFT or the c™ _1 NOT-gate 
in various coupling topologies of n qubit systems. Here 
we present examples with n up to seven. For n > 3, 
the resulting time complexities are bounded from above 
by K AK-type decompositions taken to sub-Riemannian 
regimes 

UlEES. Finally we give an outlook gener- 
alising the methods developed from spin- to pseudo-spin 
systems. 

Although the applications presented here refer to time- 
optimised quantum computing in the NMR-limit of fast 
local controls, the methods introduced are very general 
and apply to all systems whose dynamics can be cast 
into the closed form of finite-dimensional Lie algebras 
(to sufficient approximation). 



II. CONTROLLABLE SYSTEMS 
A. Spin- and Pseudo-Spin Systems 

Here we address fully controllable [U |H 03 El 
l26l ] quantum systems represented as spin- or pseudo-spin 
systems, i.e. those in which — neglecting decohcrcncc — 
for any initial state represented by its density operator 
A, the entire uni tary orbit U(A) := {UAU^ 1 \ U unitary} 
can be reached |2jfl . In systems of n qubits (e.g. spins- 



iL this is the case under the following mild conditions 
UJlliEHlh (1) the qubits have to be inequivalent i.e. 
distinguishable and selectively addressable, and (2) they 
have to be pairwise coupled (e.g. by Ising interactions), 
where the coupling topology may take the form of any 
connected graph. 

B. Time Scales for Local and Non-Local Controls 

Let the quantum system evolve in a time interval t}~ 
under combinations of piece- wise constant control Hamil- 
tonians {Hj} and the drift H^, i.e. the free-evolution 
Hamiltonian, according to 

=H A + Y,nf ) Hf ) . (1) 

3 

In NMR spin dynamics 29], for instance, the local con- 
trols on qubit t are represented by a linear combination 
of the Pauli matrices {o~i Xl o~i v }. And the drift term is 
governed by the weak scalar couplings (reminiscent of 
Ising interactions) 

H c \ = 7T Jim \ Ofe ® a mz , (2) 

Km 

provided the couplings between spins are much smaller 
than the difference between the eigenfrequcncies (shifts 
ft) of the respective spins: \Je m \ <C \tle — Q TO |. This is 
the case in heteronuclear spin systems. And in quantum 
control even the homonuclear ones can be designed such 
as to meet this greatly simplifying approximation |30| . 

For the system to be fully controllable in the sense out- 
lined above, {Hd} Li {Hj} has to form a generating set of 
the Lie algebra su(N) by way of the Lie bracket. 

Often the time scale for local controls is considerably 
faster than for the costly slow coupling evolutions. 

III. TIME-OPTIMAL QUANTUM CONTROL 

In order to control a quantum system of n qubits 
(spins- i) such as to realise a quantum gate or module 
of some quantum algorithm given by the unitary propa- 
gator Uq G U (2™) in minimal time, one has to decompose 

Uq ~ U(T) = e -^ HiM) ■ ■ ■ e-^ Hik) ■ ■ ■ e^^ (3) 

— up to a global phase factor — into a time-optimal se- 
quence (T :— J2k tk = rain) of evolutions under piece- 
wise constant Hamiltonians H ^ . 

A. Relevance of Global Phases 

However, as propagators generated by the traceless 
spin Hamiltonians are elements of the respective special 
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FIG. 1: Global phase dependence of the times needed to 
implement the 3-qubit QFT on a linear chain (L3) of nearest- 
neighbour interactions with uniform weak scalar J-couplings. 
The right curves (•) show the special unitary implementation 
of QFT with the smallest global phase <f>o — , where it takes 
2.53 J" 1 to reach a trace fidelity > 0.99999. The left curves 
(o) display the fastest QFT implementations obtained. They 
are attained with the global phase cf>i = Trace fidelities 
> 0.99999 are reached after 2.05 J -1 . Times for local controls 
are assumed to be negligible in this limit matching the typical 
NMR scenario, where the time cost is determined by coupling 
evolutions, (a) gives the trace fidelities against time, while (b) 
shows the devitiations from full fidelity in a semilogarithmic 
way. 



In Fig. this is further illustrated for the 3-qubit 
QFT implemented on a chain of three spins connected by 
nearest-neighbour interactions of weak scalar coupling in 
the NMR limit of zero time cost for the fast local controls. 



B. Optimal Control on Projective Groups 

For a given unitary quantum gate Uq and propagators 
U = U(t) describing the evolution of the quantum sys- 
tem, there are the two geometric tasks, one that explic- 
itly carries the phase, while the other one automatically 
absorbs it as desired: 

(1) minimise the distance \\U — Uq\\ 2 by maximising 
$1 := Retr{U G U}; 

(2) minimise the angle ^(U,Uq) m.od(n) by maximis- 
ing $ 2 := \tr{U G U}\ 2 . 

(1) In terms of control theory, the first task is to max- 
imise the quality functional ®i[U(t)] = Re tr{£/(l[/(T)} 
with < t < T subject to the equation of motion 
U(t) = -iHU(t) (with H = H d + Ysj UjHj) and the 
initial condition U(0) = 1, whereas the final condition 
U(T) is free at an appropriately fixed final time T (vide 
infra). As usual, the problem is readily solved by intro- 
ducing the operator- valued Lagrange multiplier X(t) sat- 
isfying X(t) = —iHX(t) and a scalar- valued Hamiltonian 
function 



unitary groups, the quantum gates Uq are realised by 
U (T) just up to global phases </> p 



U G = e-^'UiT) 



(4) 



so U G £ U(N), while U(T) e SU(N). For n spins-i read 
N := 2" henceforth. Note that with the centre of SU(N) 
being 

Z N :={e^ 1 N |p = 0,l,...,JV-l} , (5) 
one has a choice of N such phases 

P e{0 o + ^b = O,l,...,iV-i} , (6) 

where 0o shall be the smallest angle 0o G [0, ir] so that 
det{e l ^° Uq} = +1. Although global phases are clearly 
immaterial to quantum evolutions po 1— > UpoU^ 1 , it is 
important to note they do in fact contribute substantially 
to the over-all time needed to implement Uq: consider, 
e.g., 



if -if (<rfe<8>l+l®CT ms! ) 



(7) 



where the non-scalar part of the right-hand-side can be 
realised solely by (fast) local controls, whereas the left- 
hand-side hinges on nothing but (slow) coupling evolu- 
tion. 



h(U{t k )) = Re tr{At(* fc )( - i(H d + UjH s ))U(t k )} 



(8) 



Then, Pontryagin's maximum principle [31| may be ex- 
in a quantum 

dh(U(t k )) _ 



n prmcij 

ploited in a quantum setting |32l l33| to require 

-Imtr{At(t fc )fl>l7(t fc )} = (9) 



as well as the final condition for the adjoint system 

thus allowing to implement a gradient-flow based recur- 
sion |33j . For the amplitude of the jth control in iteration 
r + 1 at time interval t k one finds with a as a suitably 
chosen step size 



ff +1) (**)=«r(**) + a^ 



(r). 



.dh(U(t k )) 



(11) 



The procedure is then repeated for a set of decreasing 
final times T up to a minimal time r still allowing to get 
sufficient fidelity (compare Fig. ^| . 

(2) The second task amounts to maximising $2[U(t)] = 
\ti{U G U(T)}\ 2 , which is equivalent to the square of the 
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trace fidelity and is easy to handle by gradient-flow meth- 
ods. This problem, however, can readily be reduced to 
task (1): observe that to U G SU(N), 

U:=U*®U (12) 

is a representation of the corresponding element of the 
projective special unitary group 

SU{N) iso U{N) 



PSU(N) 



(13) 



U(l) 

embedded in SU(N 2 ). Hence this representation is 
highly reducible yet very convenient, because 

*![#(*)] = Rc tr{tf G U(T)} = |tr{^[/(T)}| 2 = <f> 2 [U(t)} . 

Hence one may adopt the previous results to obtain the 
gradient flow on PSU (N) just by using 



dh(U(t k )) 



-21m (tr{\\t k )Hft*{t k )}-tr{\(t k yU(t k )}) 

(15) 



in Eqn. 1111 Thus an explicit tensor product never enters 
the algorithm. And the final condition of the adjoint sys- 
tem does not require any prior knowledge or screening 
of the global phase ultimately giving the fastest imple- 
mentation as has been the case in previous settings, e.g. 
PI . because embedding PSU(N) in SU(N 2 ) enforces 
a global phase of zero. Absorbing the phases cuts the 
number of computations for n-qubit systems by a factor 
N = 2™. 

Having reduced task (2) to task (1) also saves all 
the convergence and step-size considerations from 
SU(N 2 ) to apply to PSU(N). 

With these stipulations, the Hamiltonians H k accord- 
ing to Eqn.^ and the numbering as in Eqn. [31 the itera- 
tions r of Ean. llll can be used in the following algorithmic 
scheme [33j : 

1. set initial controls u^\t k ) for all times t k with 
k = 1,2, . . .M at random or by guess; 

2. starting from Uo — 1, calculate for all ti, £2, ■ ■ ■ t k 
the forward-propagation 



U^{t k ) 



-i{t k -t k _ 1 )H { k r) e -i(t k - 1 -t k - 2 )H ( k r } 1 



x e -*(*i-*o)ffi r V 



(16) 



3. likewise, starting with T = and X(T) from 
Eqn. ^3 compute for all tuy ^m-i, ■ ■ - t k the back- 
propagation 

f( r ) at. . . \ wM 



x e 



A(T) 



(17) 



4. calculate dh ^ tk ^ according to Eqn. |5J 

5. with Uj +1 '(t k ) from Eqn.^Jupdate all the piece- 
wise Hamiltonians to H^ 1 ^ and return to step 2. 



Topology 
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FIG. 2: Four ordered connected graphs with 4 vertices 
representing the topology of pairwise couplings (edges) be- 
tween 4 qubits (vertices). Times given are for the shortest 
QFT-realisations obtained by numerical time-optimal control 
rounded to 0.01 J" 1 . 



IV. APPLICATIONS 

For simplicity, the coupling strengths in all the subse- 
quent examples are assumed uniform, thus enabling to 
give the times in units of J -1 . However, all our algo- 
rithms can equally well cope with non-uniform coupling 
constants directly matching the experimental settings. 



A. Towards a Time-Optimal Quantum Fourier 
Transform 

The quantum Fourier transform (QFT) is central to 
all quantum algorithms of abelian hidden subgroup type 
[35I l36| . The time required to implement this module in 
n-qubit systems clearly depends on the topology of the 
coupling interactions. 

Fig. [21 shows some of the topologies for the couplings 
of four qubits and the respective times (best numeri- 
cal results from optimal control) for implementing the 
4-qubit QFT. Clearly, the complete coupling topology 
corresponds to the maximal graph K n and thus al- 
lows for the fastest implementation, while the chain of 
nearest-neighbour interactions L n is the minimal con- 
nected graph entailing the slowest implementation. Note, 
however, that the minimal times also depend on the or- 
dering in the graph, because permutations (carried out 
by transpositions) may call for timewise costly swaps. 

The decomposition into standard gates (controlled 
phase gate, Hadamard, and swap) can readily be made 
time-optimal only in complete coupling topologies (K n ). 
There the minimal time can easily be expressed in units of 
J" 1 as a function of the number qubits (compare |37ll3^|) 



T-(n) = -(n + 3) 



(18) 



where the constant covers the final SWAP. 

However, as is shown in Fig.[3]and Tab. [I] (note the de- 
tails in the table caption), in linear spin chains (L n ) with 
nearest-neighbour Ising interactions, time-optimal con- 
trol provides a decomposition of the QFT that is much 
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number of qubits number of qubits number of qubits number of qubits 



FIG. 3: (a) Gate complexity of the QFT in linear coupling 
topologies L„. Standard-gate decomposition (•) |39ll and op- 
timised scalable gate decomposition (A) 40]. (b) Time com- 
plexity of the QFT in linear coupling topologies. Upper traces 
give analytical times associated with the decompositions of 
part (a): standard-gate decompositions (•) [39(| and opti- 
mised scalable gate decompositions (a) 40]; (a) gives a spe- 
cial (i.e. non scalable) five-qubit decomposition into standard 
gates obtained by simulated annealing Eoi . Lowest trace: 
speed-up by time-optimal control with shortest numerical re- 
alisations obtained (o) rounded to 0.01 J" 1 . Further details 
in Tab. |T| (At 7 qubits current decompositions show a trace 
fidelity of ~ 0.99 and thus have not enough significant digits 
to be included in the subsequent table.) 



TABLE I: Speed-up of the Quantum Fourier Transform on 
Linear Spin Chains, L n 



qubits stand. QFT a Blais best results c speed-ups 





T [1/J] d 


r [1/J] d 


T [l/J] d 


stand. 


Blais 


2 


1.75 


1.75 


1.25 


1.40 


1.40 


3 


8.13 


5.13 


2.05 


3.94 


2.50 


4 


17.56 


8.50 


3.15 


5.58 


2.70 


5 


30.03 


11.88(8.81) 


4.44 


6.77 


2.67(1.9 


6 


45.52 


15.25 


5.43 


8.38 


2.81 



FIG. 4: (a) Network complexity of the C n_1 NOT-gate on 
complete coupling topologies K n |4l|. (b) Time complexity 
of the C n_1 NOT-gate on complete coupling topologies. Upper 
trace: analytical times for decomposition into standard gates 
(•) H3]. Lower trace: speed-up by time-optimal control with 
shortest times (o) currently needed for realising c" _1 not by 
numerical control rounded to 0.01 J -1 . 



TABLE II: Speed-up for the C n NOT-Gate in Complete Cou- 
pling Topologies of n Qubits, K n 



qubits stand, decomposition™ best results 6 speed-up 
r [1/J] C r [1/Jf 

2 0.5 0.50 1.00 

3 3.0 1.01 2.97 

4 7.0 1.90 3.68 

5 15.0 3.37 4.45 

6 31.0 (4.59) (6.75) 



"Barenco et al. 0] 

''upper bounds to minimal time for trace fidelities > 0.99999 (for 
6 qubits currently: > 0.999) by numerical optimal control 
c times r are rounded to 0.01 J~ 1 



B. Towards a Time-Optimal c" ^OT 



" ana lytical times for decomposition into standard gates l39l 
6 |4d in brackets: the non-scalable special 5-qubit QFT 
c upper bounds to minimal time for achieving a trace fidelity of 
> 0.99999 by numerical optimal control 
"times r are rounded to 0.01 J~ 1 



faster than the corresponding decomposition into stan- 
dard gates would impose: in six qubits, for instance, the 
speed-up is more than eight-fold and in seven qubits ap- 
proximately nine-fold. 

For a fair comparison, however, note that Blais [loj 
permutes output qubits for saving SWAPS. However, 
searching through n\ permutations is beyond our cpu- 
time credits, but may provide even faster realisations in 
the future. 



Likewise, one may strive to implement the C™ _1 NOT- 
gate in a time-optimal way. In a complete coupling topol- 
ogy of n qubits, the alg orithmic complexity was described 
by Barenco et al. |4l| as increasing exponentially up to 
six qubits, whereas the increase from seven qubits on- 
wards was said to be quadratic. Again, time-optimal 
control provides a dramatic speed-up in this case as well, 
see Fig. 0] and Tab. ILTIas well as the controls in Fig. [5] 



C. Beyond Spins: Controlling Coupled Charge 
Qubits in Josephson Devices 

Obviously the optimal control methods presented thus 
far can be generalised such as to hold for systems with 
finite degrees of freedom allowing for a pseudospin for- 
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time [1/J] time [1/J] 

FIG. 5: Time course of controls for the shortest realisations 
obtained for the following quantum modules: (a) the QFT 
on a linear coupled four-qubit system (L4); (b) the C 3 NOT- 
gate on a fully coupled four-qubit system K4,. Traces in blue: 
amplitudes for cr^-controls (s-pulses); red: amplitudes for 
a"£ a -controls (y-pulses) on the spins £ — 1,2, 3, 4. 

mulation in terms of closed Lie algebras. Suffice it to 
mention the standard CNOT-gate can be realised in two 
coupled charge qubits of a solid-state Josephson device 
some five times faster than in the pioneering setting of 
Nakamura [42^. Yet one easily obtains a trace fidelity 
beyond 0.99999 as will be shown elsewhere. With the 
same fidelities one finds realisations of the TOFFOLi-gate 
in three linearly coupled charge qubits that are some nine 
times faster than by standard gate decomposition. 

V. DISCUSSION 

The goal is to extend the optimal control methods to 
larger modules of quantum algorithms or simulations in 
order to implement them both in a time-optimal and ex- 
perimentally robust way. Thus the growing set of nu- 
merical examples will hopefully provide inspiration to 
understand time-optimal steerings of quantum systems 
algebraically, which, however, seems very demanding in 
the cases presented here (compare Fig. |3J). In other in- 
stances such as for the propagator 

the theory is fully understood, and the predictions based 
on sub-Riemannian geodesies |2J| perfectly match (i) the 
time-complexity as well as (ii) the actual time course for 
the controls |£| for all irJt e [0, §]. 

Along these lines, the above controls may finally trigger 
a theoretical understanding. The ultimate challenge then 
is to extract a principle for a scalable control scheme from 
the set of numerical examples. 



VI. CONCLUSION 

Here we have left the usual approach of decomposing 
quantum modules into sets of discrete building blocks, 
such as elementary universal quantum gates thus express- 
ing the cost as algorithmic network complexity. Instead 
we proposed to refer to time complexity as the experi- 
mentally relevant cost: it allows for exploiting the con- 
tinuous differential geometry of the unitary Lie-groups 
as well as the power of quantum control for getting con- 
structive upper bounds to the time complexity both per- 
fectly matching the experimental setting while being as 
tight as possible, in particular when local and non-local 
operations are of different time scale. In the limit of zero 
cost for the fast local controls we gave decompositions for 
the QFT and the multiply-controlled NOT-gate C™~ 1 not 
that are dramatically faster than the best decompositions 
into standard gates known so far would impose. However, 
there is no guarantee the ultimate time optimum is at- 
tained, also because permutations of the qubits may give 
further improvement. 

The approach also clearly shows that in the generic 
case there is no simple one-to-one relation between time 
complexity and network complexity. This is for very 
practical reasons: typically (1) not all the elementary 
gates are of the same time cost, but each experimen- 
tal implementation comes with its characteristic ratio 
of times required for local vs non-local (coupling) op- 
erations; (2) not all the elementary gates have to be 
performed sequentially, but can be rearranged so that 
some of the commuting operations (e.g. controlled phase 
gates between several qubits) or operations in disjoint 
subspaces can be parallelised; (3) the coupling topology 
between the qubits does not have to form a complete 
graph (K„) but may be just a connected subgraph, and 
each graph comes with a specific potential of parallelis- 
ing timewise costly interactions; this is demonstrated for 
the QFT on complete coupling topologies K n versus the 
linear coupling topology L n , where the parallel perfor- 
mance of controlled phase gates [5| reduces quadratic 
time complexity to linear complexity, which, however, 
can be further speeded up by time-optimal control; (4) 
the experimental setting with its specific type and in- 
dividual strengths of coupling interaction (e.g. Ising or 
Heisenberg-Xy or XYZ type) related to the choice of 
universal gates for the network decomposition may in- 
troduce some arbitrariness. 

It is for these very reasons that time complexity is the 
more realistic measure of the experimentally relevant cost 
than network complexity is. 

VII. OUTLOOK 

Although extrapolation may be premature, it is fair 
to anticipate that in systems of some 20 qubits network 
decompositions will often become impractical. Thus 
time-optimal decompositions into controls actually avail- 
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able in the experimental setting promise to widen the 
range of experimentally accessible tasks significantly and 
will prove useful in many experimental implementations. 
Moreover, analysing the topology-dependence of minimal 
times while allowing for non-uniform coupling strengths 
will contribute valuable guidelines for designing opti- 
mised architectures of quantum computational hardware. 

By parallelising routines the results are currently being 
extended to more qubits so that time complexities can be 
deduced from fitting times against number of qubits with 
confidence. 
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